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Abstract 

We analyze the conditions for getting the Casimir repulsion between two nonequal plates. The 
■ force between plates with magnetic permeability defined by Drude or Lorentz models is calculated. 

'■^ . The short and long distance limits of the force are derived. The Casimir set-up with the hypothet- 

ical perfect matching metamaterial is discussed. We put into question the possibility of getting 
repulsion within the design of metamaterials based on metallic inclusions. 
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I. INTRODUCTION 



The study of Casimir repulsion which always presented purely academic interest has 
moved in the last years to a more practical field. It is inspired by new precision measurements 
of attractive Casimir force and the development of micro (nano) electromechanical machines 
where the Casimir repulsive force, if any, might resolve the stiction problem (see 1|] and ref. 
therein) . 

In the present paper we address the Casimir repulsion between parallel plates owing to 
the optical properties of the material. From Lifshtz formula Q it follows that the force 
between two parallel plates is repulsive if the gap between them is filled up with a material 
that satisfies the inequality £i(i£) < £3(1^) < £2(1^) [21], where £1 and £ 2 stand for the 
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dielectric permittivities of the plates, and £3 corresponds to that of the filling 
not consider this set-up in the present paper. 

The force may become repulsive if one of the plates has nontrivial magnetic permeability, 
jj, ^ 1. This possibility was not seriously regarded since for "natural" materials, where the 
magnetization of the system is due to the movement of the electrons in the atoms, yu(co>) = 1 
at visible range 5fl. However in composite materials if the inclusions are smaller than the 
wavelength, but larger than the atomic size the effective dielectric and magnetic functions 
can be introduced as a result of local field averaging. That is why the artificial materials [f| 
with magnetic response arising from micro (nano) inclusions have recently become good 
candidates for observing the Casimir repulsion. 

The Casimir repulsion for materials with dielectric permittivity and magnetic permeabil- 
ity that do not depend on the frequency was considered in [7]. The case of Drude-Lorentz 
dispersion relations for dielectric permittivity and magnetic permeability of the metamateri- 
als was analyzed in |4( . The upper limits for Casimir attractive and repulsive forces between 
multilayered structures at finite temperature were established. Starting from the Lifshitz 
formula it was shown that at zero temperature —7/8 Fc(L) < F(L) < Fc{L), where Fc is 
the force between perfect conductors. 

In jg] the Casimir repulsion due to the presence of dispersive anisotropic materials with 
gain was first predicted. These media are beyond the scope of our paper. 

In Section II we start from the basic formulas for the Casimir force between two nonequal 
plates and give the conditions for the Casimir repulsion. We consider several models for the 
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plates and get analytic results for the force at short and long distances. In section III the 
perfect matching metamaterial is discussed. We also consider recently reported optical range 
metamaterial as a candidate for Casimir repulsion. In Conclusion we discuss the possibility 
of getting repulsion within the design of metamaterials based on metallic inclusions. 



II. THE CASIMIR FORCE BETWEEN NON-EQUAL MIRRORS 

The Casimir force between two flat mirrors separated by a distance L is given by 

oo oo 

F{L) = — A Y [ duu 3 [ daa 2 - f^ B p p . (1) 

v ; 4vrc 4 ^ J J e 2aujL/c - r p A r p B K ' 

p 1 

Here r p (iu,a), p = TE,TM, are the reflection coefficients at imaginary frequencies for the 
mirrors facing vacuum 



tm V( n2 -l) + a 2 ~ea TE ^{ri 1 - 1) + a 2 - pa 

r = , r = — (2) 

a/ (n 2 — 1) + a 2 + ea y (n 2 — 1) + a 2 + pa 

with n = v /e/i, e = e(iuj) > 1, \i — p(ioo) > 1 B, ^]. The variable a is the ratio of the 
transversal wave-vector at imaginary frequency k = \Juj 2 /c 2 + k 2 to the inverse wave-length 
lo/c. 

The sign of the force is defined by the sign of the integrand in (TjQ). As |r(io;,a)| < 1, 
a "mode" {a;, a} gives a repulsive contribution to the force if the corresponding reflection 
coefficients of the mirrors A and B have opposite signs. This happens if the mirrors are 
different, ta ^ tb, and at least one mirror has nontrivial magnetic permeability. In 0, [h]] 
it was proved that no multilayered dielectric mirror can reverse the sign of the force. 

At the lower limit of the integral over a the reflection coefficients for TE and TM modes 
coincide 

lim r™ = lim r TE = {yfii- + v^)- 

a—>l q^I 

At the upper limit the reflection coefficients for TE and TM modes are different: 
lim r™ = (1 - £)/(! + e) < 0, lim r TE = -(1 - p)/(l + p) > 0. 



For fixed u one can find the values of a where the reflection coefficients change their signs: 
ocq M (uj) = y/n 2 — \ j\Je 2 — 1, ot^ E (yj) = \Jn 2 — 1/y^p 2 — 1. We split the integral over a in 
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([I]) in two: Ff(L) and F^L). In the former the integration goes from 1 to «q, and in the 
latter from a(j to oo. 

The reflection coefficients have no extremum with respect to a. If \x = 1, lim^oo r TE = 0, 
then r TE , r™ < at all frequencies. When fi < e, r™ is always negative, falling from 
r™(a = 1) to r™(a = oo), while r TE is monotonously growing from negative to positive 
values. If fi > e, r TE remains positive and growing with a, as r™ decreases from positive 
values at a = 1 to negative at a = oo. Below we analyze the force in these situations 
considering Drude or Lorentz models for the mirrors. 

A. Two non-magnetic mirrors. 

When fi = 1, neither r™ nor r TE changes the sign in the course of the integration over 
a from 1 to oo. Hence if both mirrors are non-magnetic the function ([T]) is always positive, 
and the force is attractive at all distances. As |r | > |r T£; |, the contribution of the TM 
modes to the Casimir force ([I]) exceeds the TE contribution. For the mirrors described by 
the Lorentz model Eiioj) = 1 — ul J {u 2 — Uq + i^e^u), i = A,B, at short distances L « \ ei , 
X ei = 27Tc/u ei , with uj eA < uj eB we get 



where Vt 2 u = uj 2 J2, Vt 2 2i = u 2 ei /2 + Uq } i = A,B, and G k = G k {uj Q /uj e A,ojQ/oJ e B)- The 
absorbtion in the material influences more the small frequencies which make the decisive 
contribution to the force at large plate separations. That is why the relaxation parameters 
7 ei j do not enter the short distance asymptote. 

For the mirrors described by Drude model, u Q = 0, the result is simplified to 



with G ~ 1.744, G x ~ 0.436, G 2 ~ 0.215, G 3 ~ 0.133, ... 

The long distance limit is obtained by expanding the integrand in (OQ) in powers of the 
small parameter A e ^/L or \ £ b/L and given by 




(3) 





L»X e A, A e B 



vFcas(L), r^l-4(A eA + A eB )/(37rL). 
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B. Mirror A is purely dielectric, mirror B is purely magnetic. 



Let mirror A be purely dielectric, fiA — 1- Then ,r^ M < 0. The mirror B is purely 
magnetic, Eb = 1, r^ E ,r^ M > 0. Then both TE and TM modes are repulsive at all 
frequencies, and the force is repulsive at all distances. 

The short distance limit for plasma models, Ea = 1 + ^1a/ uj2 ^ ^b = 1 +wig/w 2 , is given 

by 

The short distance attraction due to the interaction between surface plasmons is absent in 
the present case. The TM-plasmonic mode of purely dielectric mirror and TE-plasmonic 
mode of purely magnetic mirror are not coupled, and therefore do not contribute to the 
Casimir force. It results in the unusual short distance asymptote (J5]). 
At long distances we get 



L>>A e 4, A e 



vFc a s(L), ^ ~ —7/8 + 7 (\ e A + A^b) / (6irL) . 



The repulsive force coinciding with the first term of this expansion was obtained by Boyer 11] 
for two non-dispersive mirrors with ea = oo, ha = 1 and £b = 1,//b = oo. 



C. Mirror A is purely dielectric, mirror B is mainly dielectric 

Let the mirror A be purely dielectric, /j,a — 1, with r\ E ,r\ M < 0, and mirror B mainly 
dielectric, so that 1 < [ib < £b for all frequencies. Then r™ < 0, and the contribution of 
the TM modes is attractive at all distances, whereas r^ E < only for a < a$ E corresponding 
to mirror B. When a > a% the signs of the TE reflection coefficients for the mirrors A 
and B are opposite. The sign of the force is the result of the balance between F™(L) < 0, 
Ff E (L) < 0, and F? E (L) > 0. 

If Hb{u) = £b{u), then aQ E = 1, and the contribution of TE modes is entirely positive, 
F TE (L) > 0. However as r# = —r T B E < and |r^(ici;)| < |r^ E (io;)|, the total force is 
attractive (negative). Consequently, when mirror B is mainly dielectric, /xb(u;) < £b(w)> the 
force is attractive at all distances as well. At short distances it is determined by the modes 
of TM polarization and given by ([3]) . 
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FIG. 1: Reduction factor w = F/Fq as a function of dimensionless distance A = 2-7rL/A m = iv m L/c 
D. Mirror A is purely dielectric, mirror B is mainly magnetic. 

Let the mirror A be purely dielectric, fiA = 1, with r^ E , r™ < 0, and mirror B is mainly 
magnetic, so that 1 < eb < /•*£ for all frequencies. Then r^ E > 0, and the contribution of the 
TE modes is repulsive at all distances, whereas r™ > only for a < a™ ■ When a < a™ 
the signs of the TM reflection coefficients of the mirrors A and B coincide. Thus F TE (L) > 0, 
F[ M (L) > 0, F 2 ™(L) < 0. The negative term F 2 ™(L) becomes dominant for the TM modes 
at distances L < c/ {uJq a^ M ). At short distances the TM reflection coefficients are larger 
than the TE ones. The total force is attractive for short plate separation and repulsive at 
middle and long distances. For the short distance asymptote see Eq. ()3]). 

Fig. [1] gives the reduction factor of the force between purely dielectric mirror A described 
by Drude model, £a(i^) = 1 + &p/[u(l + ^)]> ^ P = 10c/L, 7 = O.Olc/L, /i^ = 1, and 
mirror B with e B (i^) = 1 + ul/(oj 2 +co% + J e &), ^fl(i^) = 1 + wi/(w 2 + uj% + 7 m u;), where 
u e = 7 e = 7 m = c/L, uo = 0.1c/ L, lc/L, 10c/ L. Varying u m L/c one can see that the 
curves cross the horizontal axis when u m ~ u e , in other words, when the material turns into 
mainly magnetic. For a metamaterial with uj m > u e , ujq/2tc ~ 3 • 10 4 GHz the force becomes 
repulsive at the distances L > 10~ 6 m. 

III. THE PERFECT MATCHING METAMATERIAL FROM THE CASIMIR 
FORCE VIEWPOINT 

If the real part of the refractive index, n = n' + in", is negative, then the transversal 
wave-vector, k z = (n 2 u 2 /c 2 — k 2 )' 1 ^ 2 , for the propagating waves k z {n',n") = —k z (—n',n"), 
and for the evanescent waves k z (n',n") = k z (—n',n"). It follows from the analysis of k z = 
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^/P\P2 exp(i ^ 1 ^ 2 ) with n' < 0. Here 

Pl = (jn'u/c- \k n \] 2 + [n"u/cfj ^ , p 2 = ([n'u/c+ \k {{ \] 2 + [n"u/cf) 1 ' 2 

n" , n" 

(pi = arctan — — — — + n, (p 2 = arctan — — — — + 

— \n'\ — \kn\- — \n'\ + Km - 

II I M I £j I I I M I ui 

When the dielectric permittivity, magnetic permeability and the refractive index tend to 
—1 [12I . [l3 |. and the absorbtion in the material is negligible, the ideal situation of perfect 
matching between material and vacuum could be achieved. It means that the transversal 
wavevector for the propagating waves k z — > —k v z ac . Consequently, the reflection coefficients 
of the propagating waves vanish, r ™^ r TE —> 0, on the interface vacuum-metamaterial. 

Let us consider two perfectly conducting mirrors one coated with a metamerial. Equa- 
tion JfBrequires that e(iu) and p{\uj) are real positive functions in accordance with causal- 
ity 5, ^]. Obviously, this condition is not satisfied for the material with constant neg- 
ative dielectric permittivity and magnetic permeability The straightforward substitution 
jj, — e — — 1 in the formula ([T|) for the Casimir force leads to positive force at L > d, and 
divergency at L < d. Moreover, the energy density of the electromagnetic field inside the 
non-dispersive material with e = fi = —1 is negative. Therefore its very existence contradicts 
the Pointing theorem, W = [d(ue(u)) /du]E 2 + [d(cop(uj))/duj]H 2 < 0. On the contrary, 
the real parts of e(u) and p{uj) may tend to —1 at a certain frequency, leaving the energy 
density of the electromagnetic field positive, W > 0. 

The Casimir energy for a multilayered system given on Fig. [2] (left) can be defined in 



terms of the scattering phase shift 5 
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f^E/^{E^,i) + /^^^). (a) 

The first term in (J6j) corresponding to bound states (surface plasmons) is absent in the case 
of perfect mirrors. 

The Maxwell equations are reduced to <f>"(z) — {/^ — EiPiOJ 2 / c 2 }(p(z) = 0, for TM modes, 
and ip"(z) — {k 2 — SiPiUJ 2 / 'c 2 }ip{z) = 0, for TE ones, % = 1..4. To get the reflection 
and transmission coefficients and the phase shift one has to solve the system of match- 
ing conditions, two for each interface in each polarization: p+ip+ = P- 1 ^--, = 
£+(fi+ = £-4>-, 4>' + = <$>'-• Here p±, e±, <fi±, ip± stand for the values of the functions when z 
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FIG. 2: The reduction factor for two mirrors, one coated by metamaterial; d = L% — L 2 



tends to the interface from the right (left). Taking the solutions for TE modes in the form 

= 1..4, B 4 — 0, ki = (ei^iuo 2 /c 2 — A^) 1 / 2 , and substituting them 



-iki 



into the matching condition we arrive at the system for the coefficients A4, Bi. Then the 
transmission and reflection coefficients are given by the ratios t(k\) = A4/A1, r(ki) = B\/A\. 
The TM transmission and reflection coefficients are obtained by replacing £j «-> //j. 

With k\ = k^ k 2 = —k 3 , e 2 = A*2 = A*i = 1, £3 = /U3 = ^3 = — 1 the scattering phase 
shift of the the four layered system is reduced to 



2i(J(fci, \a'\) = In 



t(kt 



In 



1 _ r 2 e -2ifc 2 |a'| 



(7) 



t(-ifei) ~" 1 -r 2 e 2ifc 2l«'l ' 

where a' = 2L 2 - Li - L 3 > 0, r™ = (k 2 e 1 - fci)/(fc 2 £i + A*), r TB = (fc 2 - k l )/(k 2 + fci). 
Here the limit of infinite plate separation is subtracted. If the thickness of the metamaterial 
is larger than the separation of the mirrors, a' = 2L 2 — L\ — L 3 < 0, then S(k 2 , — \a'\) = 
-6{k 2 ,\a'\). 

When layers 1 and 3 are made of perfect metal, r TE = r™ = 1, we arrive at the Casimir 
result 



,-,/ /\ , ^C7T 2 „. ,, , fiCK 2 



"720|a'| 3 ' ^ v ' 1 "~240|a'| 4 
with upper (lower) sign corresponding to positive (negative) effective distance a' and at- 
tractive (repulsive) force. The force diverges at a' = 0. This result can not be recovered as 
a limiting case of any dispersive model of a metamaterial consistent with Kramers-Kronig 
relations. 

The account for finite conductivity of the metal and dispersion in the metamaterial leads 
to a finite result at a plate separation equal to the thickness of the MM-coating. Fig. [2] (right) 
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gives the reduction factor for the force between two equal mirrors one of which is coated with 
mainly magnetic metamaterial. The mirror A and the substrate of mirror B are described 
by Drude model. The coating of mirror B is a metamaterial with effective permittivity and 
effective permeability given by Drude-Lorentz model. All parameters are normalized on the 
position of the resonance u = 27rc/A , oj p /ojq = 10, ■y/u = j e /u) = Im/^o = 0.01, 
uo e /ujQ = 2, uj m /ujQ = 3. Here d = L3 — L2 is the thickness of the metamaterial. For 
comparison we plot the reduction factor rf = F(a')/Fc(a'), a' = (L/\q — 1) for two perfect 
metals one coated by perfect matching metamaterial. 

We see that if the mirror B is two-layered, the force changes the sign two times. For 
short and intermediate distances, L < d, the contribution of the reflections from the in- 
terface between the layer of the metamaterial and the substrate is small, and the behavior 
corresponding to bulk MM-mirror is reproduced (compare with Fig. [T|) . At large distances, 
L > d, the fluctuations "feel" the presence of the substrate, and the reflection coefficients 
approach the ones for the substrate. The force becomes attractive, achieving the values 
typical for metals. When the thickness of the metamaterial is smaller than its characteristic 
wavelength Ao, the layer becomes transparent for fluctuations with uo > ujq, and the region 
of repulsion is considerably narrowed (solid curve in Fig. [2]). 

Further we discuss recently reported low loss optical metamaterial For describing the 
material we use the effective media approach, considering anisotropic compound material 
as a homogeneous media having effective dielectric and magnetic functions. To evaluate 
the appropriate parameters of the effective media we took the complex permittivity and 
permeability plots from Fig. 3 in 161 ] . 

For the dielectric permittivity we have taken the Drude response with a small anti- 
resonance, the magnetic permeability has Lorentz form, Fig. [31 



uj {uj + ifl e ) uj 2 — cUq + iVL e iU ' uj 2 — oj 2 + ifimicu ' 

with the position of the resonance defined by ujq. Though the amplitude of the magnetic 
resonance is larger than the corresponding dielectric one, the Drude-type term dominates 
in e(u). The magnetic permeability of the meta-material is described by a function which 
is characteristic for the dielectric permittivities of the semiconductors, for example sili- 
con (FigJHD- 

The calculation of the Casimir force between golden mirror A and mirror B made of 



9 




FIG. 3: Left panel: e(iw) for meta-material discussed in [16| in comparison with gold and silicium. 
Right panel: the magnetic permeability of the meta-material fJ,(iu); fJ,Au = ^Si = 1- 

this metamaterial was performed with the following values of the parameters: io e j2 r n = 
3.6 ■ 10 5 GHz, Q e /2n = 8.9 ■ 10 3 GHz, u /2tt = 2.05 • 10 5 GHz, u el /2n = 2.04 • 10 4 GHz, 
fi e i/2vr = 5.03 • 10 3 GHz, u ml /2n = 9.72 • 10 4 GHz, Q ml /2n = 1.1 • 10 4 GHz. The 
metamaterial being mainly dielectric, we predict no repulsion in the present setup. 

The effective medium approach is valid for wavelengths longer than the "lattice constant" 
of the meta-material. To put it differently, the theoretical estimations for the force are 
trustable for plate separations large in comparison with the "lattice constant" of the meta- 
material. For more accurate results optical data in a wide frequency range and for different 
incidence angles are needed. 



IV. CONCLUSION 



When one of the mirrors is mainly magnetic and its magnetic permeability is described 
by the Drude-Lorenz model the repulsion could be observed at the separations of the order 
L ~ Ao = 2-kc/ujq. However the force decreases rapidly with the distance. That is why to 
get measurable Casimir repulsion we look for a material with magnetic permeability /i ^ 1 
at optical frequencies. 

Non- magnetic conductive elements, such as split ring resonators (SRR), being embedded 



into a dielectric lead to nontrivial magnetic response of the compound [17| . A metamaterial 
with effective permeability /i e // is obtained when a lot of such elements are oriented in 
different directions and positioned as cubic lattice [6]. As the effective resonance frequency 
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is defined by the geometric scale s of the resonator, loq ~ 1/yLC ~ 1/s [18|], the latter 
should be of several hundred nm size. Though when the frequency of the field approaches the 
plasma frequency of the metal this estimation is not valid anymore, as the electrons become 
insensitive to the variation of the field. Experimentally it was shown that geometrical scaling 
law for the resonance frequency brakes down at about 300 THz [19I ]. 

We seek for mainly magnetic material. The height of the magnetic resonance peak for a 
single resonator increases with the filling fraction. On the contrary, the interaction between 
the inclusions makes the resonance peak broader and lower. At the same time the metallic 
inclusions change the dielectric permittivity of the host media. It acquires the properties of 
diluted metal or highly doped dielectric with the dielectric permittivity dominated by the 
Drude term. In other words, providing us with needed magnetic response the metallic struc- 
tures add as well redundant dielectric permittivity, that makes the metamaterials mainly 
dielectric. Consequently they do not manifest repulsion in the Casimir set-up. The meta- 
materials based entirely on dielectrics are more appropriate. That could be, for example, 
arrays of dielectric spheres in a dielectric matrix [20j 



, |21( with the magnetic response arising 



from polaritonic resonances. 
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